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Here a viable and never before investigated mechanism of propagation of spikes along neurons 
is proposed. In the following, velocity of propagation of small-amplitude pressure waves through 
the cytoplasmic interior of myelinated and unmyelinated axons of different diameters is theoretically 
estimated, and is found to generally agree with the action potential (AP) conduction velocities. This 
remarkable coincidence allows to surmise a model in which AP spread along axon is propelled not 
by straggling ionic currents as in the widely accepted local circuit theory, but by mechanoactivation 
of the membrane ion channels by a traveling pressure pulse. Hydraulic pulses propagating in the 
viscous axoplasm are calculated to decay over ~ 1 mm distances, and it is further hypothesized 
that it is the role of influxing during the AP Ca 2+ ions to activate the membrane skeletal protein 
network for a brief radial contraction amplifying the pressure pulse and preventing its decay. 

The model correctly predicts that the AP conduction velocity should vary as the one-half power 
of axon diameter for large unmyelinated axons, and as the first power of the diameter for myelinated 
axons, provided that specific mechanical properties of axons are independent from diameter; that 
myelinization increases the conduction velocity; that the conduction velocity increases with the 
temperature. Unlike the local circuit theory, the model is able to qualitatively explain observed 
increase in the AP duration in axons of smaller diameters. Predictions of absolute AP conduction 
velocities are limited by the knowledge of relevant to propagation of pressure waves mechanical 
properties of axons, still, the velocities are predicted well for myelinated axons, while an agreement 
for unmyelinated axons requires 3 orders of magnitude higher resistance of axons to deformation 
increasing their diameter, compared to values deduced from published data on membrane area 
expansion moduli. Experimental test of the model is needed. 



Depolarization of interior of an excitable cell to a 
critical level leads to spiking behavior of voltage across 
the cell membrane, or action potential (AP). The spike 
has been shown to arise from rapid changes in mem- 
brane ion-specific permeability, allowing flow of trans- 
membrane ionic currents powered by electrochemical gra- 
dients 0. Membrane permeability to ions was found 
to be regulated by voltage across membrane, and a self- 
regenerating scheme of time course of the AP was de- 
veloped 0- The rapid changes in membrane perme- 
ability were later shown to arise from brief activation 
of voltage-gated ion-selective protein channels embedded 
in the membrane These findings are supported by 
more than half a century of experiments and are a basis 
of modern research in neuroscience Q . 

Phenomenon of propagation of the AP through ex- 
citable cells (as opposed to development of the AP at 
a given cell point) has been explained by the local circuit 
theory, which postulates that spread of ionic currents 
powered by AP-associated voltage spike across mem- 
brane depolarizes adjoining unexcited membrane and 
brings it to the critical level for excitation 0|. Predic- 
tion of velocity of AP conduction is an important test of 
this model, and can be found from solution of a second- 
order partial differential equation for time-dependent po- 
tential along the cell length, with parameters such as 
electrical permeability of cytoplasm and membrane, and 
membrane capacitance 2J. A simplified form of this 



'Electronic address: maratr@mit.edu, rvachev@alum.mit.edu 



equation was numerically solved in J2| for a squid giant 
axon model, with membrane permeability to Na + and 
K + parametrized as a function of voltage history across 
membrane, as deduced from voltage-clamp data. The 
conduction velocity was found to be 18.8 m/s for a 476 
/im diameter axon at 18.5 °C with 21.2 m/s experimental 
value. 

The local circuit theory predicts that conduction ve- 
locity in unmyelinated nerve fibers should vary as the 
one-half power of axon diameter if axoplasmic resistiv- 
ity and membrane electrical properties per unit area are 
constant, with experimental data fitting powers between 
0.57 and 1 depending on the fiber class [lj. For myeli- 
nated fibers, the local circuit theory predicts that the 
velocity should vary as the first power of diameter, again 
with provisions of constant specific electrical properties, 
and with additional constraint of proportional scaling of 
axon diameter, external myelin diameter, distance be- 
tween the nodes of Ranvier and membrane area at the 
nodes 4]. The first power relationship is well established 
for large myelinated fibers 0. 

Experiments show that in smaller diameter axons of a 
related class, duration of the AP spike becomes longer, 
a fact completely unexplained by the local circuit theory 
. Conduction velocity increases with the temperature, 
e.g. for cat vagus myelinated fibers Qio (ratio of the 
velocity at one temperature to the velocity at a tempera- 
ture 10 °C lower) is 4.8 at 18 °C, 2.5 at 28 °C and 1.6 at 
37 °C and for desheathed rabbit vagus unmyelinated 
fibers Qi is 3.5 at 10 °C, 2.1 at 20 °C and 1.7 at 30 
°C 0. These changes are explained by the local circuit 
theory only qualitatively, through increased rate of con- 
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formational changes in the membrane channel proteins 
at higher temperatures 

Overall, except for the numerical calculation of the 
conduction velocity for the squid giant axon model 0, 
agreement between local circuit theory predictions and 
measurements is close, but rather loose. Agreement in 
dependence of the velocity on axon diameter for unmyeli- 
nated fibers requires variation in specific electrical prop- 
erties of axoplasm or membrane with diameter, while for 
similar agreement for myelinated fibers these are assumed 
to be constant with (experimentally not strictly satis- 
fied) assumptions of geometrical scalings [lj. Tempera- 
ture affects conduction velocity similarly in myelinated 
and unmyelinated fibers, while if due only to changes in 
activation rate of membrane channels, simple reasoning 
suggests that myelinated fibers would be affected much 
less, since their conduction velocity would change at the 
Ranvier nodes only. As already mentioned, longer du- 
ration of the AP in fibers of smaller diameters is not 
explained. 

The purpose of this paper is to show that by making 
two assumptions, namely, that a hydraulic pulse prop- 
agating through axoplasm along axon length is able to 
mechanoactivate enough membrane ion channels for the 
influxing currents to depolarize membrane to the excita- 
tion level, and that at least one of influxing during the 
AP ion species, presumably Ca 2+ , activates the mem- 
brane skeletal protein network for a brief contraction 
timed to amplify the propagating hydraulic pulse, one 
obtains a competing with the local circuit theory model 
of AP propagation, and that this model agrees reason- 
ably well with experimental data on the AP conduction 
velocities. 

The exact mechanisms of how a traveling hydraulic 
pulse may be amplified and how it might elicit membrane 
depolarization are uncertain and most appealing to au- 
thor processes are discussed here. However, exact nature 
of these processes does not change basic model features 
and conclusions drawn, as long as these processes do am- 
plify the hydraulic pulse and do excite the membrane 
upon arrival of the pulse. Just as in the local circuit the- 
ory, it is assumed that depolarization of a membrane seg- 
ment to the critical for excitation level causes membrane 
voltage to spike due to activation of membrane voltage- 
gated ion channels. The distinction from the local circuit 
theory is that this initial depolarization is brought about 
by a propagating hydraulic pulse and not by the strag- 
gling ionic currents, and thus velocity of AP conduction 
is equal to velocity of propagation of the hydraulic pulse. 
It should be noted that correct prediction of AP shape by 
the Hodgkin-Huxley equations does not immediately 
refute the present model, since local development of AP 
is the same in both mechanisms. 

Although most of discussion here is concerned with the 
AP propagation in neurons, the model is easily general- 
ized to non-neuronal excitable cells. 



I. MODEL FEATURES AND NECESSARY 
ASSUMPTIONS 

Hydraulic pulses do propagate through viscous liquids 
enclosed in flexible tubes, a fact illustrated by blood pulse 
propagating through blood vessels. In addition to nor- 
mal force on the tube wall exerted by the pulse, viscosity 
of the liquid introduces shear stress in the direction par- 
allel to the axis of the tube [|| . These forces acting on an 
extensible axolemma should lead to a membrane stretch, 
which, according to the first assumption, alone or in com- 
bination with bending and compression of the membrane 
by the propagating pulse, should open enough ion chan- 
nels to bring membrane voltage to the excitation level. 
Evidence is now mounting that many unrelated types 
of ion channels are (for unclear reasons) mechanosensi- 
tive 0, but still it is not here possible to quantitatively 
substantiate this assumption, partly because magnitude 
of the membrane deformations depends on an unknown 
pulse amplitude. Since many mechanosensory ion chan- 
nels are activated by protein links anchoring channel pro- 
teins to extracellular structures and to the cytoskeleton 
|8|, it is also reasonable to propose that perturbation 
of the cell interior cytoskeleton by the propagating hy- 
draulic pulse could be directly mechanically coupled by 
the protein links to the channel proteins, for a fast and 
sensitive activation. 

As will be shown below, pressure waves in axoplasm de- 
crease e-fold in amplitude over ~ I mm distances, so that 
their sustained propagation in long axons would require 
amplification. It appears that forceful circumferential 
(and consequently radial) contraction of the membrane 
skeletal protein network, triggered by influxing during 
the AP Ca 2+ ions, is the most plausible mechanism of 
such amplification. Indeed, virtually every excitable cell 
has voltage-gated Ca 2+ channels in its membrane and, in 
fact, calcium is the basis of the AP in muscles of many 
invertebrates, in smooth muscles and in many gland cells 
of vertebrates |3|]. At the same time, increase in free in- 
tracellular calcium is often associated with initiation of 
motion, from motility in freely moving cells and mus- 
cle contraction to synaptic vesicle release in synapses . 
Influxing during the AP calcium would instantly many- 
fold increase its concentration in the membrane skeletal 
protein network near the inner membrane surface and 
that, coupled with calcium ability to quickly induce con- 
formational changes in proteins, such as upon its bind- 
ing to actin filaments in actin-myosin muscle complex, 
could provide for a fast contractile response necessary to 
amplify a quickly propagating pressure pulse. The well- 
known sodium impulse in this scheme could simply be the 
means to rapidly change transmembrane voltage by the 
sodium current and thus quickly open relatively sparse 
voltage-gated calcium channels. 

The purpose of the voltage spike in this model is to 
activate voltage-gated calcium channels, but clearly the 
spike will create axoplasmic ionic currents that, just as 
in the local circuit theory, will depolarize adjacent unex- 
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cited membrane. This depolarization, then, will decrease 
number of mechanoactivated by the hydraulic pulse ion 
channels needed to bring the membrane to the excitation 
voltage; but if the straggling currents are able to excite 
adjacent membrane before arrival of the pressure pulse, 
AP will obviously spread according to the local circuit 
theory mechanism. Thus, it appears that the two mech- 
anisms of propagation can be viewed as competing, and 
it is conceivable that different cells, depending on their 
particular electrical and mechanical properties, are able 
to realize one or the other mechanism. 

Another argument in favor of pressure pulse ampli- 
fication by structures adjacent to the membrane inner 
surface is the existence of (currently functionally unex- 
plained) dense "undercoating" just beneath axon mem- 
brane in the axon hillock and in the nodes of Ranvier [j| . 
Initiation of a hydraulic pulse in the axon hillock, as well 
as restoration of its amplitude in the nodes of Ranvier af- 
ter passive propagation of the pulse through myelinated 
part of axon, would require extra effort, and it is rea- 
sonable to expect that features underlying pulse creation 
and amplification should be emphasized. Expression in 
these regions of the dense membrane undercoating might 
then suggest that structures adjacent to the membrane 
inner surface take part in creation and amplification of 
the pulse, possibly, as hypothesized above, by a forceful 
circumferential contraction of a protein network included 
in the undercoating. Another well-known feature of the 
nodes of Ranvier, high density of voltage-gated sodium 
channels, can be justified within the model by necessity 
to rapidly change transmembrane voltage by the sodium 
current to quickly open voltage-gated calcium channels, 
by an increase in membrane depolarization at the next 
Ranvier node due to augmented axoplasmic currents, and 
possibly by an increase in number of sodium channels 
mechanically coupled to interior cytoskeleton structures, 
thus improving mcchanosensitivity of the membrane to 
axoplasmic disturbances. 

In the model being presented axon acts as a sen- 
sor of mechanical disturbances in the axoplasm, trans- 
ducing them into pressure-generating motion by utiliz- 
ing mechanosensitivity of the membrane ion channels 
and ionic electrochemical gradients across the membrane. 
Considering preponderance of mechanical stimuli in na- 
ture, an appealing scheme of how cells could have evolved 
to transmit signals by pressure waves can be constructed. 

Let's suppose that a cell has developed an ability to 
rapidly withdraw a part of its membrane to which an 
external negative pressure and/or stretch was applied. 
This reaction could be aimed at detachment from ad- 
hesion surfaces, or simply at preservation of shape of a 
freely moving cell in a dynamic mechanical environment. 
If such cell happens to have cylindrical shape, and an 
external stimulus causes rapid membrane withdrawal at 
a cell end, a wave of increased cytoplasmic pressure will 
be created and will propagate along the axis of the cell 
toward its other end. Since passive membrane movement 
is determined by pressure difference across it, circular 



membrane segments subjected to increased intracellular 
pressure in the pulse will move as if negative pressure 
was applied at the extracellular side; membrane stretch 
caused by the motion of the viscous cytoplasm in the 
pulse should also be indistinguishable from a membrane 
stretch exerted from outside the cell. If these stimuli are 
sufficient to elicit forceful inward retraction of the circu- 
lar membrane segments, new cytoplasmic pressure waves 
will be created, which, if timed properly, will amplify the 
original hydraulic pulse carrying the signal of stimulation 
at one cell end to the other. 



II. PASSIVE PROPAGATION OF HYDRAULIC 
PULSES THROUGH AXONS 

For the case of small-amplitude harmonic waves in in- 
compressible liquid enclosed in a thin- walled elastic tube, 
with the wave wavelength large compared to the radius 
of the tube, and the case of "large" liquid viscosity (de- 
fined below) , speed of propagation of the waves Vi nc and 
decay length Li„ c (length over which amplitude of the 
waves decreases e-fold) are |(J 
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where E is the Young's modulus of the tube wall material, 
h is the thickness of the tube wall, R is the radius of the 
tube, p is the liquid density, /i is the liquid viscosity, uj is 
the wave frequency, v is the Poisson's ratio of the tube 
wall material. 

The condition of "large" viscosity (or small tube ra- 
dius, or low frequency) is R(ujp/ p) 1 ^ 2 <C 1. That this 
condition holds well for axons of diameters at least up to 
~ 40 pm can be seen from the following estimate. Mea- 
surements of macroscopic cytoplasm viscosity with ~ 1 
/im objects freely diffusing or moving under ~ 1 pN forces 
yield /i ~ 0.2 Pa s [rHHI, while application of larger ~ 
1 nN forces to objects of similar sizes yields p ~ 210 
Pa s 01 • Duration of the pressure pulse should be close 
to that of the contraction that created and amplified it, 
which in turn should be determined by duration of influx 
of calcium ions, rate of dissipation of free intracellular 
calcium and time properties of the presumed process of 
contraction. Lacking detailed knowledge of the last two 
processes, it is here simply assumed that duration of the 
pressure pulse is equal to the AP duration, e.g. ~ 0.34 
ms for large-diameter cat myelinated axons at 37.1 °C 
P3, with corresponding uj w 2tt/(2 • 0.34 ms) » 9200 
rad/s. Taking R = 20 /Ltm, p = 1000 kg/m 3 , p = 0.2 
Pa s, uj = 9200 rad/s, yields R{ujp/pf/ 2 w 0.14, which 
should be smaller enough than 1 for the large viscosity 
limit to hold. Another condition used to derive (1) - (2), 
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of large wavelength A of the pressure waves compared to 
the radius of the tube, or A = 2irv/u ^> R, as will be 
seen later, also holds well for axons. 

For the case of small-amplitude harmonic waves in 
a compressible liquid enclosed in a rigid tube, and the 
above conditions of large viscosity (R(upj p) 1 ! 2 -C 1) and 
large wavelength (A > R) 



in a compressible liquid enclosed in a thin-walled elas- 
tic tube, with the conditions of large liquid viscosity and 
large wavelength: 
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where v soun d is the speed of sound in the liquid in the 
absence of viscosity. 

It should be noted that in the absence of viscosity, 
speed of propagation of small-amplitude harmonic pres- 
sure waves in a compressible liquid enclosed in a thin- 
walled elastic tube is given by (e.g. 0) 
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where k is the intrinsic bulk compressibility of the liquid, 
and k + 2R/Eh is the compressibility of the liquid in the 
tube stemming from both the intrinsic compressibility of 
the liquid and distensibility of the walls. In the limit of 
rigid tube walls so that their distensibility can be ignored, 
i.e. when k > 2R/Eh, (5) reduces to w M=0 = (1/pk) 1 / 2 , 
which is equal to the speed of sound in the liquid v soun d 
(e.g. 0) and the first factor in (3) and the second in 
(4). In the opposing limit of k <C 2R/Eh, i.e. when 
internal compressibility of the liquid can be ignored, (5) 
simplifies to v M =o = (£'ft./2i?p) 1 / 2 , which is known as the 
Mocns-Korteweg equation and the first factor in (1) 
and the second in (2). It now can be seen that outcome 
of introduction of large liquid viscosity in the inviscid 
case is independent from whether compressibility of the 
liquid in the tube arises from distensibility of the walls, 
or from intrinsic compressibility of the liquid, except for 
an additional factor 2/(5 — 4^) 1 / 2 in (1) and (2), which 
takes into account movement of the tube wall in the di- 
rection of the tube axis by the viscous liquid . In prin- 
ciple it is a straightforward, although cumbersome, task 
to modify equation of conservation of mass (3) in [(| to 
include intrinsic compressibility of the liquid in the ap- 
proximate form dp = kpdp, repeat derivations and obtain 
expressions for the wave speed and decay length for the 
case of a compressible viscous liquid in an elastic tube. 
Here, however, it is noted that the factor 2/(5 — Av) 1 / 2 
varies from 0.89 to 1.15 for physical values of Poisson's 
ratio < v < 0.5, and is close enough to 1 to be ne- 
glected in the present analysis. Therefore, neglecting ax- 
ial motion of the tube wall, (1) - (5) can be combined to 
obtain for the case of small-amplitude harmonic waves 
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In general, when velocity of harmonic waves (phase 
velocity) v depends on frequency of the waves u>, velocity 
of propagation of a wave packet is equal to the group 
velocity v gr , given by (e.g. 
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Application of (8) to (6) yields group velocity v gr = 1.5?;. 

Derivations in Q assumed isotropy of the tube wall 
material, while structure of lipid membranes suggests 
high anisotropy of their elastic properties in the trans- 
verse direction compared to the in-plane directions. By 
following the derivations, however, one can ascertain that 
terms corresponding to the transverse elasticity of the 
wall material have to be neglected. Then Young's mod- 
ulus E and Poisson's ratio v in the formulas above rep- 
resent membrane resistance to a tangential to the mem- 
brane stretch and shrinkage of perpendicular to stretch 
in-plane dimensions respectively, under assumption of in- 
plane isotropy of membrane elastic properties. Then, ne- 
glecting transverse stresses on the membrane, Hook's law 
yields (e.g. [13): e\ = (o\ - va 2 )/E, e 2 = (cr 2 - vo\)jE, 
where e\, e 2 - relative orthogonal elongation of mem- 
brane subjected to orthogonal tangential stresses cri, er 2 . 
From the definition of elastic area expansion modulus K 
(e.g. [3): (ti+t 2 )/2 = KdA/A, where Ti and r 2 are or- 
thogonal tensions in the membrane surface, dA/A - frac- 
tional area change. Taking into account that n = <J\h, 
t 2 = a 2 h, and for small ei, e 2 : dA/A = e\ + e 2 , one can 
obtain: 



Eh = 2K(1 - v) 



(9) 



III. NUMERICAL ESTIMATES AND 
COMPARISON OF MODEL PREDICTIONS TO 
EXPERIMENT 

Inspection of formulas (6) - (8) in relation to propa- 
gation of pressure pulses through axons leads to the fol- 
lowing statements: velocity of propagation v for a given 
u, /i, k, R increases with increased membrane rigidity 
Eh; for a given uj, p, k, Eh, dependence on the axon 
radius is intermediate between R 1 ^ 2 and R and specifi- 
cally depends on the product Ehk; for a limiting case of 
a soft membrane (Eh <C 2R/k), v ~ R 1 ^ 2 ; for the op- 
posing case of a rigid membrane (Eh 3> 2R/k), v ~ R; 
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for an axon of fixed R, Eh and k: v ~ (w/^) 1 / 2 ; decay 
length L behaves similarly to v, except for the L ~ w -1 / 2 
dependence. 

As has been mentioned earlier, the local circuit the- 
ory is unable to explain increased duration of the AP in 
the smaller diameter fibers of a related class |l|. Since 
in the model being presented membrane channels are ac- 
tivated not only by voltage across membrane (as in the 
local circuit theory), but also by the traveling pressure 
pulse, hydraulic pulses of longer duration should increase 
duration of the AP due to longer mechanostimulation of 
the membrane channels. Since lower frequency compo- 
nents of a pressure pulse decay less with distance (from 
(7)), increased damping in fibers of smaller diameters 
(again from (7)) in general means that propagating pres- 
sure pulses in these fibers will have lower frequency and 
longer duration, and, as mentioned above, longer dura- 
tion of the AP. 

From the above considerations it also follows that du- 
ration of a "membrane" AP in this model is predicted 
to be independent from the fiber diameter, and equal to 
values obtained with the local circuit theory, owing to 
absence of a propagating pressure pulse in this case. Du- 
ration of a propagated AP, however, is predicted to be 
longer than that for a membrane AP, and is predicted 
to increase with decreasing fiber diameter reflecting in- 
creased duration of the pressure pulse. The idea that 
the AP duration ATap adjusts to the duration of the 
pressure pulse AT pr reinforces the approximate relation 
AT^p ~ AT pr , assumed earlier from converse arguments 
that duration of the pressure pulse should be close to 
that of influx of calcium ions. 



A. Myelinated axons 

As follows from (6) - (8), the larger the stiffness of 
the axolemma, the larger the speed of propagation of 
the pressure waves and the lesser their decay with dis- 
tance. From this perspective, myelinization of axons 
by Schwann or oligodendroglial cells tightly wrapping 
around axon in a spiral manner is aimed at increase of 
rigidity of the axon membrane. Successive layers of the 
same myelin cell are known to firmly adhere together 
by proteins from the cell surface glycoprotein family, 
which arc in fact main protein constituents of myelin 
0. This adherence should be effective against unwind- 
ing of the spiral when subjected to increased pressure 
in the hydraulic pulse, and should transfer the distend- 
ing load to lipids and proteins of all myelin layers. In 
order to be able to neglect distensibility of the myelin 
sheath compared to intrinsic compressibility of the ax- 
oplasm, "effective" Young's modulus E of the myelin 
should satisfy stemming from (6) approximate condition 
E > 2R/(hk). Taking k as for saline water at 37 °C, 
k = {pv 2 ^^)- 1 « 4.04-l(r 10 Pa" 1 0, and R « h leads 
to the condition E > 5 • 10 9 Pa. Area expansion modulus 
K for plasma membranes is on the order of 0.4 N/m with 



membrane thickness of h w 6 nm [18J , which corresponds 
to Eh s» 0.6 N/m (from (9), v taken as 1/4) and "effec- 
tive" Young's modulus of E = Eh/h « 10 8 Pa. Thus, if 
myelin sheath is as elastic as plasma membrane, it can- 
not be considered indistensible compared to the intrinsic 
compressibility of axoplasm. However, considering that 
lipid composition of myelin differs markedly from that of 
plasma membranes, in particular in that it contains sig- 
nificantl y m ore cholesterol known to increase membrane 
rigidity |l8t l20|. that the protein-rich major dense lines 
might also significantly contribute to myelin rigidity, and 
that in the model being presented increased myelin stiff- 
ness is beneficial for propagation of the pressure pulses 
and thus should be sought after by the cellular mecha- 
nisms, it is here assumed that circumferential distensibil- 
ity of myelin sheath does not appreciably effect velocity 
of propagation of axoplasmic pressure waves. Then, also 
neglecting reduction of the propagation velocity in the 
nodes of Ranvier, the quotient 2R/(Eh) in (6) and (7) 
can be disregarded for myelinated axons. 

Substituting in (6) compressibility of saline water k = 
4.04 • 10~ 10 Pa -1 , viscosity fx — 0.2 Pa s, frequency 
u> = 27r/(2 • 0.6 ms) » 5200 rad/s, axon radius R = 1 
[Am leads to the pressure pulse velocity v gr = 6 m/s, 
which, given the approximate nature of the estimate, is 
in a good agreement with the experimental value of 10 
m/s for the AP conduction velocity in cat myelinated 
axons of 1 /im radius at 38 °C jj; the AP duration of 
0.6 ms for the mentioned fibers was measured at 37.1 
°C in 13] . Use of fx = 0.2 Pa s in the estimate can 
be justified by considering that it should represent vis- 
cosity arising from small-amplitude cytoskeletal deforma- 
tions 0, E| that are likely to accompany propagation 
of a small-amplitude pressure pulse, while larger values 
of fi = 210 Pa s were obtained with forceful displace- 
ment of cytoplasmic objects that possibl y d isrupted the 
cytoskeleton or distorted it considerably [lj] ; still, given 
the large scatter of values for ll available in the literature, 
uncertainty in ll should represent a major contribution 
to overall error of the estimate. 

From (7), e-fold amplitude decay length for the con- 
sidered case is then L w 1.2 mm, and wavelength A = 
lixvjio Ri 7.3 mm. 



B. Unmyelinated axons 

For unmyelinated axons of radius R = 0.65 /im, as- 
suming membrane Eh = 0.6 N/m as before, and taking 
fx = 0.2 Pa s, k = 4.04 • lO" 10 Pa" 1 , u = 5200 rad/s as 
for the myelinated axons, (6) and (8) yield the pressure 
pulse velocity v gr = 0.053 m/s, which is ~40 times lower 
compared to measurements of 2.3 m/s for the AP con- 
duction velocity in cat unmyelinated nerve fibers of 0.65 
Lira radius at 38 °C Q], notwithstanding the low value for 
viscosity used. Taking into account agreement obtained 
with the same values of fx, k, ui for myelinated axons, 
and that the main difference in mechanical properties 
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between myelinated and unmyelinated axons should be 
in the axolemmal resistance Eh to increases in axon di- 
ameter, one can conclude that in order for the model to 
predict correct AP conduction velocity for the unmyeli- 
nated axons, "effective" membrane Eh has to be ~ 1600 
times larger than the value of 0.6 N/m used. This lat- 
ter value represents resistance of lipid plasma membranes 
to area expansion. Additional resistance to increases in 
diameter of unmyelinated axons can come from defor- 
mation of adjacent glial cells; from deformation of the 
cytoskeleton tethering the membrane; from passive re- 
sistance of the membrane skeletal protein network to the 
radial deformation and from active contraction of the net- 
work during amplification of the hydraulic pulse. Since 
axial membrane elasticity does not appreciably influence 
velocity and decay length of the pressure waves, the men- 
tioned above factors might selectively increase membrane 
resistance to the circumferential expansion, leaving ax- 
ial rigidity at a lower value. From this perspective, the 
model predicts "effective" area expansion modulus for 
the deformation described ~ 3 orders of magnitude larger 
than for lipid membranes. 



C. Velocity dependence on axon diameter and 
temperature 

From (6) - (8) it follows that for unmyelinated fibers of 
small diameters (2R <§; Ehk), v gr ~ R, while in the op- 
posing limit of large diameters (2R 3> Ehk), v gr ~ i? 1 / 2 . 
In unmyelinated C-fibers (diameters 0.4 - 1.2 fim) of cat 
saphenous nerve AP conduction velocity vap is propor- 
tional to R, in unmyelinated fibers (1.6 - 20 fim) of locust 
and cockroach vap ~ R ' 7 - R 08 , while in unmyelinated 
fibers (2 - 520 fim) of cephaloid molluscs vap ~ R 057 
0. Thus, assuming that Ehk is on the order of 10 jum 
and is not significantly different for the mentioned fibers, 
predicted gradual change in dependence from R to R 
with increasing unmyelinated fiber diameter is very well 
matched by experiment. 

For myelinated fibers, assuming Eh ^> 2R/k, (6) - 
(8) give v gr ~ R dependence, again closely matched by 
experimentally determined vap ~ R Q- In the con- 
sidered above limit of small unmyelinated fibers, condi- 
tion 2R -C Ehk is the same as for myelinated fibers, 
Eh ^> 2R/k, which indicates that for these fibers myelin- 
ization would not increase the speed of propagation or 
decay length. Very small diameter fibers are indeed usu- 
ally unmyelinated Q. 

The analyses above implicitly assumed that specific 
mechanical properties of axons, k, fi, Eh, are constant 
across fibers of varying diameters, or that their varia- 
tion is small enough not to appreciably change obtained 
R dependences. The same condition was assumed for 
the central frequency lj of the pressure wave packet. 
In cat saphenous myelinated fibers at 37.1 °C, AP du- 
ration ATap changes from 0.6 ms to 0.34 ms with 8- 
fold increase in diameter |l3[ . Corresponding change in 



lj = 2tt/(2 • AT A p), from 5.2 • 10 3 to 9.2 • 10 3 rad/s, can 
be shown to introduce not more than an additional 0.14 
power of R in the predicted dependences. A larger vari- 
ation of lu with R, however, would have a greater impact 
on the predictions. 

Assuming that Eh, k and R do not significantly depend 
on temperature T, and that changes in to and fi with T 
are independent from degree of fiber myelinization, (6) 
- (8) predict that myelinated and unmyelinated fibers 
should be similarly affected by a temperature change, 
just as observed experimentally [l|; the predicted depen- 
dence is vap ~ (w(T)//i(T)) 1 / 2 . Experimental data on 
fi(T) is not readily available, however, for most liquids 
viscosity increases with decreasing temperature [2lJ • An 
increase in viscosity means (see (7)) augmented damping 
of pressure waves of all frequencies, and, just as in the 
case of fibers of smaller diameters considered earlier, du- 
ration of the AP can be argued to increase. AP duration 
indeed increases with decreasing temperature, e.g. with 
Qio of 3.4 at 37 °C for cat vagus myelinated fibers [l3| . 
but this large variation could mainly stem from change 
in rate of activation of voltage-gated membrane channels 
PJ. According to v gr ~ w 1 / 2 and the assumed relation 
lj = 2n/(2-AT A p), Qio of 3.4 for the AP duration ATap 
implies that Qio for the pressure pulse propagation ve- 
locity is predicted to be 3.4 1 / 2 w 1.8 without even consid- 
ering fi(T) dependence. Experimental value of Qio = 1.6 
for the AP conduction velocity in the cat vagus myeli- 
nated fibers at 37 °C thus might indicate that viscos- 
ity of the axoplasm does not change appreciably in the 
27 - 37 °C range. 

It is also interesting to note that cooling-related re- 
moval of AP conduction block in damaged (e.g. by mul- 
tiple sclerosis) myelinated axons p2| can be naturally ex- 
plained within the model by increase in the decay length 
L with the decrease of lu ~ (ATap)^ 1 at lower tempera- 
tures. 



D. Depletion of sodium and calcium concentrations 

Depletion of sodium concentration around a squid gi- 
ant axon leads to wider and slower propagating APs |23j , 
in close agreement with the Hodgkin-Huxley equations 
[24|. However, longer duration of the AP in the pre- 
sented model also decreases velocity of AP propagation, 
through relation vap <~ (ATap) -1 / 2 ; qualitative analysis 
of vap and ATap from Fig. 4 in [23 does not exclude 
validity of this relation. 

Conduction in axons in low external calcium is poor 
or blocked at turnon of cathodal current |lj. These ax- 
ons, however, do respond better to anodal break exci- 
tation 1], which does provide higher gradient for Ca 2+ 
to enter the axon. Of interest is whether axons prop- 
agate APs in extremely depleted external calcium. If 
they do, interesting is to know whether conduction ve- 
locity increases or decreases, the latter being predicted 
by the model, if conduction can switch to lower velocity 
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local circuit mechanism. In the local circuit theory, con- 
duction velocity with depleted external calcium should, 
in the first approximation, increase: decrease in free ex- 
tracellular calcium significantly decreases thresholds of 
activation of both sodium and potassium channels |24j : 
free intracellular calcium can activate polarizing calcium- 
dependent potassium currents [24[. Then, observed in- 
crease in the velocity would favor the local circuit the- 
ory, but, strictly speaking, would not disprove the model, 
since it is possible that the velocity given by the local 
circuit theory became larger than that predicted by the 
model after modification of channel gating, and the lo- 
cal circuit theory mechanism took over the propagation 
from the pressure wave. Data on the AP conduction 
velocities in axons with manipulated calcium concentra- 
tions is not readily available, since most experiments use 
voltage-clamp methods under these conditions; detailed 
analysis of these phenomena will be done in a subsequent 
work. 



IV. DISCUSSION 

It is very remarkable that transmission of hydraulic 
pulses through axoplasm, a process so disparate from 
passive spread of ionic currents, reasonably well repro- 
duces all major features of AP conduction: changes in 
the velocity with membrane rigidity, with axon diameter, 
with the temperature and with AP duration; changes in 
the AP duration with fiber diameter. Accurate predic- 
tions of absolute AP conduction velocities by the model 
are limited by the knowledge of relevant to propagation 
of the pressure pulses cytoplasm viscosity, resistance of 
myelinated and unmyelinated axons to increases in di- 
ameter, cytoplasm compressibility, and by approximate 
nature of the estimate of duration of the hydraulic pulse. 
Still, estimates of the velocity for myelinated fibers, with 
an assumption of indistensibility of the myelin sheath, 
give values close to those observed experimentally; for 
unmyelinated fibers, although the velocity is underpre- 
dicted, the situation might change as the appropriate 
experimental data on membrane distensibility becomes 
available. 

It also has been noted (Prof. M. Kardar, 2003, private 
communication) that the lipid membrane can be more 
realistically modeled as 2-dimensional liquid rather than 
elastic solid, with wave propagation velocity shaped by 
bending membrane rigidity rather than area expansion 
modulus. With the different treatment, which will be 
attempted in the future, predicted by the model AP ve- 
locity in unmyelinated axons might change dramatically. 

On the other hand it is noteworthy that all fundamen- 
tal ingredients necessary for realization of the described 
model are in place: mechanosensitivity of many mem- 
brane ion channels (e.g. 0), low free calcium concen- 
tration in a resting cell, ubiquity of the calcium impulse 
and its trait to induce motion in filaments 3], and high 
filament content of the cell interior yl . Whether or not 



these ingredients are present to a degree sufficient for the 
active propagation of the hydraulic pulses, and whether 
any cells in fact use the described mechanism for AP 
conduction, should be best settled by experiment. 

A small axon swelling followed by a smaller shrinkage, 
both roughly synchronous with propagating APs, was 
demonstrated in crab and squid giant axons [25j . Beyond 
increased pressure inside axon, the swelling can be caused 
by variation of the voltage across membrane through the 
ensuing change in difference in surface tensions between 
the two membrane interfaces |26| . Axial motion of an 
axoplasmic object, such as a vesicle, observed simulta- 
neously with the AP, should be less influenced by the 
voltage-dependent membrane movement and more sensi- 
tive to axoplasm disturbances induced by the hydraulic 
pulse. The model predicts that a swelling should precede 
the voltage spike, and, in fact, figures in [25| indicate that 
the swelling by several milliseconds does precede the volt- 
age spike, but the authors' conclusion was that cause of 
that might be in electronics delays; careful analysis of 
the delays in an experiment similar to that in |25| could 
provide another test of the model. 

Knockout of a neurofilament subunit gene has been re- 
ported to ~ twice decrease conduction velocity in mice 
large myelinated axons, with axon diameter, AP ampli- 
tude, duration and shape preserved [2tJ . Neurofilament 
proteins are the most abundant cytoskeletal element in 
large myelinated axons [23] • These data are favorable for 
the presented model, with the changes in conduction ve- 
locity explained by modification of cytoskeletal and mem- 
brane skeletal protein mechanics. The local circuit the- 
ory predicts that conduction velocity is independent from 
mechanical properties of axons. Somewhat contrary to 
that, it is known that cell pressurization does block AP 
conduction in at least in some axons 0. 

Other tests can come from measurements of velocity of 
propagation of pressure pulses through axons; from pre- 
cise measurements of mechanical properties of neurons 
and other excitable cells; from experiments on mechan- 
ical stimulation of AP in excitable cells; from measure- 
ments of the conduction velocity with depletion of free 
intracellular calcium, or with disruption of the membrane 
skeletal proteins; or from a setup aimed at direct detec- 
tion of pressure waves accompanying APs. 

The following circumstances can be recognized to in- 
directly support the model: in many non-neuronal ex- 
citable cells AP seems to be connected to mechanical mo- 
tion 0] ; neurons often have convex shape as if pressurized 
from inside; myelinization appears to be the only feasi- 
ble way to substantially increase mechanical resistance of 
axolemma, while resistance of the membrane to leak cur- 
rents could in principle be increased by other means, e.g. 
decrease in number of leak channels; many freely moving 
cells have an ability to control membrane movement by 
forces of microfilaments located just beneath the plasma 
membrane 0; the same ability is attributed to neurons, 
where the forces are believed to arise from interactions 
of actin, fodrin and other elements 0. Finally, the mam 
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components of the model, propagating hydraulic pulses 
and radially contracting tubes generating a directed ax- 
ial force on the tube contents, are exploited in cardiovas- 
cular, intestinal and other tubular systems in many or- 
ganisms, which might even suggest an evolutionary link 
between the systems and excitable cells. 

Central role of hydraulic pulses and pressure- 
generating processes in the model implies that pressure- 
related phenomena should play important roles in the 
function of those neurons and excitable cells, that would 
realize the model mechanisms: synaptic vesicle release in 
neurons might be aided by impact of the AP-associated 
hydraulic pulse on the vesicles or by contraction of the 
bell-shaped membrane skeletal protein network of the 
synapse, creating a directed push on the cytoplasm and 
the vesicles; depolarization of the axon hillock to the ex- 
citation level might be accomplished not only by pas- 
sive spread of ionic currents from postsynaptic sites, but 
also by mechanoactivation of the membrane ion chan- 
nels by cytoplasm displacement or pressurization; there- 
fore, synaptic integration might involve relay of displace- 
ment or pressurization of cytoplasm through the den- 
dritic tree to the axon hillock, and learning and mem- 
ory in neurons might be partially reflected in strength 
and duration of Ca 2+ -activated contraction of postsynap- 
tic membrane skeletal proteins, or in another pressure- 
generating process; AP generated in a neuron might in- 



fluence processes in the cell body, dendrites and con- 
tacting synapses by mechanical disturbances introduced 
by the back-propagating pressure pulse; AP-associated 
hydraulic pulses might participate in transport of sub- 
stances in excitable cells. Other outcomes of the model 
are the hypotheses on the evolution of the signal trans- 
mission in cells, on reasons for mechanosensitivity of 
many unrelated types of ion channels, and on the physi- 
ological role of axonal calcium influx. 

In the case the presented model is not experimentally 
confirmed, calculations in the paper show that velocity 
of propagation of ~ 1 kHz pressure waves through axons 
can be close to the AP conduction velocities, and hence 
pressure waves generated by voltage-induced membrane 
movement during AP conduction |26| might accumulate 
into a larger shock-like wave. 
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